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Abstract 

Elementary features of galileon models are discussed at an introductory level. Following a simple 
example, a general formalism leading to a hierarchy of field equations and Lagrangians is developed for 
flat spacetimes. Legendre duality is discussed. Implicit and explicit solutions are then constructed and 
analyzed in some detail. Galileon shock fronts are conjectured to exist. Finally, some interesting general 
relativistic effects are studied for galileons coupled minimally to gravity. Spherically symmetric galileon 
and metric solutions with naked curvature singularities are obtained and are shown to be separated from 
solutions which exhibit event horizons by a critical curve in the space of boundary data. 
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1 Introduction 



Galileon theories are a class of models for hypothetical scalar fields whose Lagrangians involve multilinears 
of first and second derivatives, but whose nonlinear field equations are still only second order. They may 
be important for the description of large-scale features in astrophysics as well as for elementary particle 
theory [12] US]. Hierarchies of galileon Lagrangians were discussed mathematically for flat spacetime in 
[T51 [TOl I2"U] , independently of an earlier systematic survey of second-order scalar-tensor field equations in 
curved 4D spacetime [23]. The simplest example involves a single scalar field, 0. This galileon field 
may be coupled "universally" to the trace of the energy-momentum tensor, 0, and upon so doing, it is 
gravitation-like by virtue of the similarity between this universal coupling and that of the metric to 9 
in general relativity. As might be expected from this similarity and the ubiquitous generation of scalar fields 
by the process of dimensional reduction, it is possible to obtain some galileon models from limits of higher 
dimensional gravitation theories. Indeed, galileon models were discovered yet again by this process [13] . 



2 A simple example 

Although higher derivative actions usually lead to higher derivative field equations, as is well-known, nonethe- 
less it is possible to accommodate higher derivatives in the action while retaining second-order field equations 
if the Lagrangian is not quadratic in the fields. The price to be paid is that the second-order field equations 
are always nonlinear. This is the basic ingredient that underlies all galileon models. 

This point has already been well-appreciated in the literature, of course, but for purposes of illustration, 
consider the "simplest" cubic, fourth-order Lagrangian density, 

L = 4> 2 Qa /3/3 , (1) 

where Q — d<j)/dx a , etc., and repeated indices are summed using either Lorentzian or Euclidean signatures. 
The corresponding action A = J L has local variation, 

8 A 

— = 20 QQ(9(3 + (0 2 ) qq/3/3 . (2) 

Surface terms, irrelevant for the field equations in the bulk, have been discarded. That is to say, 

8 A 

— = 40 (j) a al3f3 + 84>a 0Q/3/3 + 20 aa 4>t3p + ^4>u(i 4>aj3 , (3) 

and the field equation, 8A/5<f> — 0, is both nonlinear and fourth-order. 

However, the order of the field equation may be reduced by adding to the Lagrangian a judicious amount 
of the "next-to-simplest" cubic term that involves first, second, or third derivatives. So far as the field 
equations are concerned, there is actually only one other term that can be added, namely, cj> aa fipp- 
Superficially different terms, e.g. <Pa <Pa/3i3, 4> 4>aj3 "Pap, 4>a 4>a 0/30, and <p a 4>p <fi a p, do not give independent 
contributions to the local variation of the action in the bulk, although they differ in their surface contributions. 
In particular, 

2 4> 4>aa 0/3/3 + 40 (j) a i3 Q/3 - 3 2 aa /3/3 = ( 4 </> 0q </> q/ 3 + 2^,3 Qa - 3 (f? cj) aa p - 2 (j) a (j) a (j>p) ^ , (4) 

i.e. a total divergence. Thus it is sufficient to include in the action any two of the three terms on the LHS, 
with an arbitrary relative coefficient. 

So, rather than ([1]), consider instead the Lagrangian density 

L = <f> 2 (f> a aPP - A (j) (f>aa 4>PP ) (5) 

with constant A. The variation of the action obtained from ([5]) is 

^ = 20 aQ/3 £ + (0 2 ) QQ/3/3 - A ^0 QQ ((>0I3 + 2 (0 0aa)^) 

= (4 - 2A) (0 <j)aaPl3 + 20 a (frapp) + (2 - 3A) Qa <j)pp + 40a/3 4>a0 ■ (6) 

Thus A = 2 uniquely eliminates from the variation all derivatives higher than the second, leaving just 

SA 

= -4 (0 QQ 0/3/3 - Q /3 4>afs) ■ (7) 

A=2 
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While this equation is still nonlinear, it is now only second-order. 

Moreover, the action for the A = 2 model can be rewritten in various ways upon integrating by parts. 
Perhaps the most compact and memorable of these is 



A 2 = j 4> a (f) a (j>pp d n x . (8) 
This differs from the previous ^4| A=2 by a factor of 2 and a boundary term, 

M = ~ J {4> 2 4>euxpp - 2 4> 4> aa <f>pp) d n x - i J d a B a d n x , 

involving the current 

B a — cf> 2 C>2 0/30 = 4> 2 <i>app -2 (j) (j) a (j)pp . (9) 

Indeed, most discussions of this model are developed around A 2 , after defining the system's Lagrangian to 
be 

L% = (j) a 4> a 4>pP • (10) 

To complete our discussion of this elementary case, consider the energy-momentum density arising from 
L 2 . The canonical result is straightforwardly obtained, even though the Lagrangian involves higher deriva- 
tives, but the resulting density is not a symmetric tensor. However, minimal coupling to gravity is guaranteed 
to yield a symmetric tensor, so we take that route. Generally covariant forms of (jTUJ) and © are obtained 
through the replacements (f) a <t>a. d n x — > g a ^ *<j> a <t>pyf—g d n x and <f)pp — > d fl (y rr gg^<j)v)- Thus an 

invariant action is 

' g a ?<p a <f>p {V^gg^fa) d n x . (11) 



^ I curved space 

Varying the metric gives Q a p. In the flat-space limit, the result is 

GV = <t>ti<t>v<t>aa - fia^aufifi ~ <t>a<l>an<t>v + &v.v4>a4>P§<xP ■ (12) 

This is seen to be conserved 

d^O^ =^£2, (13) 

upon using the field equation, £2 = 0, where 

£2 = 4>aa4>Pp - 4>ap4>al3 ■ (14) 

The justification for the name "galileon" is now apparent. Any shift of the field by a constant, or by a term 
linear in x, as is reminiscent of a galilean transformation in classical mechanics, will leave the action (J5J 
invariant, up to surface terms, and therefore not falsify a solution of the field equation. 

An interesting wrinkle now appears: M „ is not traceless on-shell. Consequently, the usual form of 
the scale current, x a Q atl , is not conserved. On the other hand, the action © is homogeneous in <j> and 
its derivatives, and is clearly invariant under the scale transformations x — > sx and <j)(x) — » s^ 4 ~ n ^ 3 (j) (sx). 
Hence the corresponding Noether current must be conserved. This current is easily found, at least in four 
dimensions, since the trace is obviously a total divergence in that case: 

®wln=4 = 4>aa4>p4>P + ^4>a4>p4>aP = 9 a {(j)a<t>p<t>p) ■ (15) 

That is to say, for n = 4 the virial is the trilinear V a = 4>a.<t>p<j>p- So a conserved scale current is given by 
the combination, 

^ln=4 = ^060^ - (j) a 4>a4>pL ■ (16) 

However, the virial here is not a divergence modulo a conserved current. So the theory is not conformally 
invariant despite being scale invariant. 

Some additional algebra is needed for n^4, but eventually one finds: 

Sfx = x a Q afl Vfi , (17) 

Tl — 1 4 — Tl 

= — — <$>ol4>cl$h h — — <i>J^ , (is) 



n — 4 

Xa 0q H o 



£2 ■ (19) 
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The last term in the virial for n^4 involves a bilinear current which is conserved merely as a restatement 
of the field equation: 

Jfi = (frfiOa&a , £2 = <9 M J p . (20) 
In fact, this current is itself a total divergence, 

Jp = % (Snu4>(t>aa ~ 4> 4>iiv) , (21) 

so the field equation for the model is a double divergence for any n. But once again, although the model is 
scale invariant in any number of dimensions, it is not conformally invariant for n > 2. 

Since (fTU)) has the form of the conventional free field Lagrangian density x the Klein-Gordon equation, 
for a massless scalar field, it immediately suggests a generalization to a hierarchy of such systems, where the 
Lagrangian density for the kth system is just a product of the free field Lagrangian density and the equation 
of motion for the (k — l)st system. In fact, this simple generalization is easy to formulate in explicit detail. 
A systematic theory for the hierarchy is elegantly expressed using determinants. 



3 General formalism 

This section may be skipped by anyone with a phobia for determinants, and definitely should be passed over 
by anyone under a doctor's orders to cut back on tensor index shuffling. Later sections of the paper rarely 
invoke results obtained in this section. However, the material presented here may be helpful for applications 
beyond those considered in this primer. Accordingly, the last part (§3.9) provides an encapsulation of the 
results in terms of determinants and standard Kronecker symbols. 

3.1 Determinant and trace identities 

For any n x n matrix M, consider the expansion 



det (1 + AM) = £k (M) , 



(22) 



where £q = 1. Elementary cases are £\ = Tr (M) and £ n = nl det (AI). Other cases may not be so familiar. 
However, from the identity 

det (1 + AM) = exp (Tr In (1 + AM)) , (23) 
it follows that the £k obey a recursion relation for any M, 



' 1 (—l) k ~ 1 ~ e 



1=0 



Tk-i£i , for k < n , where T m = Tr (M m ) . 



The solution of this recursion for all k < n can be expressed in terms of another set of determinants, 

£ k = det (T fe ) , 

where Tfc is an auxiliary k x k matrix containing the various traces: 



(24) 



(25) 



T fe = det 



/ 71 

r 3 

Tk-2 

Tfc-i 

V % 



k-1 
7! 
T 2 

Tk-3 
Tk-2 

Tk-i 





k-2 
Ti 

Tk-i 

Tk-3 
Tk-2 



\ 










Ti 2 

r 2 7i 1 
r 3 r 2 t ) 



(26) 



The recursion relation (|24p is recovered by expanding det (Tfc) in the minors of the first column. In addition 
to (|25p we also note the identity 

fef fe _ 1 = Tr(adjT fe ) . (27) 
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In this last trace relation, we use the adjugate (a.k.a. the classical adjoint) matrix notation, adj (T) 
(detT) T" 1 . 



For example, 

£i = dot Ti = 71 , 



det To = det 



T?-T 2 



(28) 
(29) 



71 1 

r 2 Ti 

71 2 

det T 3 = dot ( T 2 71 1 ) = Ti - 371 7~ 2 + 2T 3 , 

r 3 r 2 71 

/ 71 3 \ 

r 2 ti 2 

r 3 r 2 ti 1 

V T 4 r 3 r 2 Ti / 

etc. Actually, for an n x n matrix M, explicit computation of the traces inserted into the expression (1261) 

gives a tim'Z determinant for k > n. That is, £k> n = 0, as would be expected from the expansion of 
det (1 + AM). 

Moreover, a slight modification of the auxiliary matrix in (|26j) gives directly the characteristic polynomial 
for any n x n matrix M , 



dot T4 = det 



71 4 - 67?T 2 + 871T 3 + 3T 2 2 - 671 



(30) 



(31) 



/ 1 

A 
A 2 



det (M - A 1) 



det 



n 
71 

72 





n — 1 
71 



\ 





71 2 

r 2 ti 1 

r 3 r 2 Ti / 



(32) 



A" -2 7^-2 71,- 3 
A' 1 1 71,-i 71i_ 2 
\ A 71 71i 71,-1 

This follows immediately from expanding in the minors of the first column, using (|22[) and ([24] 
From the elementary identity det (AB) = (det ^4) (det B) we have 

det (1 + AM) = A" (det M) det (l + A^M" 1 ) , 

for an n x n nonsingular M . It follows for < k < n that 

^ £ k (M) = (detM) , 1 £„_ fe (Af- 1 ) . 
k\ (n — k)\ 

For example, 

£ n {M) =n!(detM) , 
£ n _i (M) = (n - 1)! (detM) Tr (M -1 ) = (n - 1)! £1 (adjM) , 

etc. Or, to rewrite (|34j) more symmetrically, for n x n nonsingular M, 

1 1 



1 i £ fc (M) = 
^det (M) fe! 



^/det (M- 1 ) (n-fe)! 



(33) 
(34) 



(35) 
(36) 



(37) 



3.2 Field equations 

Take M to be H = dd(f>, the Hessian matrix of second partial derivatives of cj> 1 

det (1 + \dd</>) = TT £ * ( dd ^) ■ 



, x„), then 



fe=0 



(38) 



For k > 1 also define "the equation of motion at level k" as £/. (dd<j)) = 0. For example, with <f) a p — d a dp<t>, 
£ x {dd<t>) = Tr (990) = c/> aa , (39) 
£ 2 {dd<t>) = (Tr {dd^)f - Tr ((dc^) 2 ) = - <f> a p<f> a p , (40) 
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etc., while at the highest levels, in n dimensions, 

£ n {dd(f>) = nl det (dd<j>) , (41) 
£ n -i (dd<p) = (n - 1)1 det {dd<j>) Tr ((dd^ 1 ) ={n- 1)! £ x (adj (dd<p)) , (42) 

etc. We shall refer to the £ n (dd(f>) — case as the "maximal" galileon field equations. 
3.3 Lagrangians 

These may be defined recursively and yield the above equations of motion after varying <j) and integrating 
by parts: 

C k = ct> a c}> a £ k -i {dd<P) , sjc k d n x = -2 J £ k {ddcfy Sep d n x . (43) 
Thus, stationarity of the action for C k implies the equation of motion: 

= £ k (dd^) . (44) 
A systematic method to obtain this recursion is to work out the variation of 

A = J <j) a (f> a det (1 + Xddcp) d n x , (45) 

and then use (f22j) to single out the action for C k . More generally, with 

A [F] = J F {(f> a (f> a ) det (1 + \ddcf>) d n x , (46) 



we find 



5 A [F] = -2 J £ [F] Sep d n x , (47) 



where 

£[F]=det(l + Xdd4>) [(<t> aa -\{l + \dd<t>)-l <P^ a )F' {1 + Xddcj>)-1 d„(4>M)F"} . 

(48) 

Setting F (<p a 4> a ) = (f'a't'a, i.e. F' = 1 and F" = 0, and expanding the RHS of (|4"5)l in powers of A leads to 
©. 

3.4 Universal field equations 

Given in n dimensions an arbitrary Lagrangian dependent only upon first derivatives of the field, 0, and 
homogeneous of weight one, there is an iterative procedure for calculating a sequence of equations of motion 
which always terminates with the same final equation, lA n = 0, independent of the starting Lagrangian 
[2171 [T7] . This final equation has therefore been called a "universal field equation" (UFE) It involves only 
first and second derivatives of 4>. Here we describe the relation between U n and the galileon Lagrangian C n 
in n dimensions. 

The functional form appearing in the UFE in n dimensions can be expressed as a "bordered determinant" 

B3I3], 

U n [4>] = det ( ° ) , (49) 



d(/> dd(j> 

where the entries in the top row, and in the left column, are and <p a for a — 1, • • • ,n, and the Hessian 
matrix occupies the n x n block on the lower right. However, unlike £ k (dd<p) that appears in the galileon 
field equations, U n is not identical to a total divergence, so the integral J U n [0] d n x can serve to specify 
nontrivial dynamics in the bulk. 



1 In two dimensions, the UFE is just the Bateman equation, <j>xx<t>t ~ '^4>xt<t>x<Pt + 4>tt<t>t = 0- 
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For example, in n — 2 Euclidean dimensions, 



U 2 




(50) 



But then a a /30£ = dp (0a0a0,3) - 0a/30a0/3 - 0a 
U 2 = -|0a 0a 0/3/3 + ! #/3 (0a 0a 0/3)- That is to say. 



022 - 01202) + 02 (02101 ~ 0201l) 

SO 0a0a/30/3 = \dp (0a0a0/3) ~ 50a0a0/3/3- ThuS 

(51) 



3 1 



The same result applies in spaces with Lorentz signature, when repeated indices are summed with the Lorentz 
metric. 

Similarly, in n = 3 dimensions, U3 differs from a constant multiple of £3 just by a divergence, 



2 ^7 (0aa 







(52) 



Indeed, it turns out that in any n dimensions, U n is always proportional to the maximal galileon Lagrangian 
C n modulo a divergence, or boundary term, i3„H 



(n - 1)! U n = ~ (n + 1) £„ + \(n- 1) B„ 



(53) 



The relative coefficient between 1l n and C n is worked out explicitly in the next subsection, where an explicit 
form for B n is also given. The upshot is that the action for maximal galileon fields that vanish on the 
spacetime boundary is obtained just by integrating the functional form appearing in the UFE, J U n [0] d n x. 

As someone well-schooled in determinants might guess, especially in light of the discussion following 
there is another way to express the UFE in n dimensions, in terms of traces. This is given by 



V„ = det 



/ So 


n-1 





■ 





\ 


Si 


Ti 


n-2 ■ 


• 








S 2 


To 


Ti ■ 


• 








S n -3 


T n -3 


Tn-4 


■ Ti 


2 





S n -2 


Tn-2 


7~n-3 


■ r 2 


7T 


1 


\ S n -i 


Tn-l 


Tn-2 


■ % 


T 


Ti J 



= - (n - 1)! U n 



(54) 



where we have defined 



T k = Tr 



{dd(t>) k 



S k = Tr 



(0090) (dd(t>y 



(55) 



The special results in ([51]) and (|52j) may be confirmed from (|54|) . 

Perhaps the most elegant proof that the determinants in (j49|) and (|54|) are proportional is to make use of 
an orthogonal transformation at each point x to find local frames such that the symmetric Hessian matrix 
is diagonal, 

(000)=diag(A 1 ,--- ,A„) . (56) 

(We do not diagonalize the full (n + 1) x [n + 1) matrix appearing in (|4l?|) because we wish to keep track of 
the first derivatives, 00.) In such frames it is straightforward to show, from either (T4T))) or (f5"4"]l . that 



U n = - 



E^ 2 

a=l 



a (11^) ■ 



(57) 



/3=1 
/3^a 



det 



Note this would still be true if the zero in the upper left corner of lA n were replaced by any constant c, for then 
( c dd 



d(j> dd<t> 



■ Un + -T- En and the last term is again a total divergence. 
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The signs here are for n Euclidean dimensions. For spaces with Lorentz signature, the result holds upon 
making the usual sign changes. 

The determinant in (I5~il) should be compared to that in (|25p . This suggests that we consider more 
generally, for < k < n, the determinant 



V k = det 



Expand this determinant in the minors of the first column to obtain 



/ S 


k - 1 





■ 





\ 


Si 


Ti 


fc-2 • 


• 








s 2 


T 


Tx ■ 


• 








Sk-3 


Tk-3 


Tk-i 


• T 


2 





Sk-2 


Tk-2 


Tk-3 


■ T 


7] 


1 


V Sk-i 


Tk-i 


Tk-2 


■ % 


T 


Tx I 



(58) 



V k = S £k-i - (k - 1) SxEk-2 + {k — l)(k — 2) S 2 £ k -3 -+••• + (-if 1 (* - l)!S*-i£b ■ (59) 

This in turn suggests a more direct — but somewhat brute force — derivation of the relation between (p[9"| 
and (|54)l . and the result in (|53| . This is given in the next section. 



3.5 More on determinant identities 

We have defined £k (M) by the expansion (|22|) . On the other hand, we also have the identities 



det (1 + AM) 
1 



I °Qi ■ ■ -a 
Oil " ' a 



m 
1 



: : xnAM Qlft+ ^:" 



X 



n(n— 1) 



A 2 Af Ql( 9 1 M Q2( g 2 + 



?QlQ2 v n(ra-l) \2 

A 1, 
n! 



»! + (n - l)!^ 1 x n AM aift + (n - 2)!d££ x 2i!fiZ A 2 M aift M Q2& + ■ 



= 1 + A S;i x M Qlft + - Sftg x M aiPl M a2h + ■ ■ ■ + — ^.VX" x M aif3l ■ ■ ■ M anPn 
where the generalized Kronecker symbols are defined bjj^| 



m r .. ttt _ I tti-a t q w -a, 



[n-k)\ 
1 

(n-k)\ 



(60) 



(61) 



= SplSp* ■ ■ ■ 5p* *8p* =t permutations of as or /3s, but not both. 

For emphasis, on the RHS we have underlined the repeated indices that are implicitly summed. Thus the 
Ek (M) can be expressed in terms of generalized Kronecker symbols: 



£ k (m) = ifzi: >< m ^ ■ ■ ■ M akPk . 

Note the last term in the expansion ([60|) is the familiar 

£ n (M) = ^X-X" x M aifil ■ ■ ■ M an f) n = n\ det M . 

3 There is a difference here between Euclidean and Minkoski metrics. For Euclidean space, 
is true for any n, but the corresponding identity in Minkowski space is 
The remaining discussion in this subsection will be given for the Euclidean case. 



(62) 
(63) 
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Applying generalized Kronecker symbol methods to bordered determinants, and making use of (|()2"j) . leads 
to relations that may be usefully applied to the UFE. Consider 



M = 



v 
v M 



(64) 



where M is any symmetric n x n matrix, v is an arbitrary n x 1 column matrix, and its transpose v is a 
1 x n row matrix. Clearly det M is bilinear in the components of v. For convenience, we index the rows 
and columns of M. from to n. Since Afoo = it follows that 

detM= J^TTy. 6 m> :: ^ x Ma °*> ' 



i 



(65) 



where other subscripts can not appear because of the antisymmetry of the Kronecker delta. So, substituting 
v and M for the components of M, 

-1 



det M = TT7 <5 fl 1 V"-/3rT-ik n x "ftV^mft ' ' ' Af a „_ l/3 „_ 1 (66) 



(n 


-1)! 




-1 


(n 


-1)! 




-1 


(n 


-D! 



detTW = 



Repeating the steps of the last three lines to reduce the remaining Kronecker symbol gives the series 

-1 [ (vv) £ n -i (M) - (n - 1) (vMv) £„_ 2 (M) + (n - 1) (n - 2) (?M 2 u) £„_ 2 (M) 
(n-1)! [ + (-l)"" 1 (n-l)!(wAf"- 1 «)£ (M) 

— 1 11-1 C _ 1 M 

TV E TrV^ (^) £ n-i-j (M) 



(67) 



where (dMS) = Tr (wM fe ) = w a (M fe ) a/3 u^, and of course, £ (M) = 1. 

In particular, setting v — d(j> and M = <93</>, and recalling the second definition in ([5"5j) . this last result 
leads directly from (|4"9")| to 



Un 



(n - 1)! ^ 

y ' 3=0 



(68) 



But the sum on the right-hand side here is just another way to write (|58|) for k — n, as expressed in (|59[) . 
Thus we establish again the relation between U n and V n , ([54^1 , only this time without invoking a local frame 
to diagonalize the Hessian matrix. 

Next, we consider the local variation of J U n [(/)] d n x using (|66|) . Thus 



U n W\ d n x = - — J Splp^J™_*p x <j>p<j) a x 0„ 1(8l • • • ctn _ 1J g n _ 1 d n a: , 
and dropping surface terms in the variation, 
8 f U n [<p\ d n x 

S Ti02---l3n~ip X ( 2 (<%) ^« X "?Wi ' ' ' 0q„_i/3„_! + (n - 1) </>/3</> Q X (5(f) aih ) (j) a ^. 



(69) 





-1 


(w 


-1)! 




-1 


(w 


-1)! 




-1 


(n 


-1)! 


n 


+ 1 



rOltt 2 -0„. 



(n- 1)! 



2 ' ' ' <Pa„_!/3„_ 

X ( - 2 ^a/3 x ' ' ' 0a„_i/3„_i + (n - 1) {4>l3<t>a) a ^ X <?W 2 ' ' ' <?W_i/3„_ 

X ^^■■■Pn-iP X ( _2< ^a/3 X ^aift ' ' ' 0an-i,9»-l + (« - 1) ^ai/^afr x ^asft ' ' ' 0a„_i/3„-i 
X ^Ifa'-'Pn-iP X ( ^"1/31^2,32 •• •0a„_i/3„_i0a„ / 3„ ) d™Z (70) 
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Now use £ n (ddcf>) as given by to obtain the result, 

8 Ju n [<t>] d n x = " + * J £ n {ddcfy 5<f> d n x . (71) 
Compare this to 8 J C n d n x = —2 J £ n (dd(j>) 8(f) d n x, as given in (|43|) for k — n, to conclude 

<5 1 Z4 [0] = 8 J (~ 2 " n + _\y £n [0]) rf"^ , (72) 

where spacetime boundary terms have been dropped. Thus the unvaried integrands can only differ by a 
divergence. This establishes (|53|) . An explicit form for the divergence B n , as normalized in (|53[) . can be 
found by keeping track of the discarded boundary terms produced by integrating by parts in (|43[) and in 
(|70| . This is left as an exercise. The result is 

&n = 9 a ( Sppipl'.'.'p^ X </>aiPi ■ ■ ■ </>a„-2/3„-2 X 0p0 7 7 ) ■ (73) 

For example, this reduces to the divergence terms in (|51|) and ([52]) . for n — 2 and n = 3, respectively. 

To complete our discussion of the UFE, we consider the effects of a quadratic constraint on cf> a for a 
field in n + 1 dimensions. The constraint will effectively reduce the number of dimensions to be n. For 
convenience, we let indices a, /3 = 0, 1, • • • , n, while we let A,/x, v = 1, • • ■ , n. Then the constraint of interest 
to us is 

<j> a <p a = . (74) 

Nontrivial solutions would of course require complex fields in the Euclidean case, but for the time being, let 
us not be deterred by this.. Solve for and differentiate to obtain 

0o = -0m0m I 00m = 7 — ) wo = -2 ■ (75) 

00 00 

Now compute £k {dd(f>) subject to the constraint, specifically displaying the occurrences of </>om and </>oo- 
Thus 

£fc (990)| n+ l dimensions = ^02-'.pk X 0<*101 ' ' ' 0afc/3fc 
with 00,00—0 

= k8 v l V2 ... v *_i X 0^!!/! ' ' ' 0M(s_lI/fc_l x 000 

- A; (& - 1) 8^'.::^ x • • ■ M(J _ 2 „ fe _ 2 x O a0oa 

' "l/jUS — V^Ml^l Yfik^k 

= k8 VlV2 ... Vk _ 1 X 4 > a 1 u 1 ' ' ' 0Mfc-l"fc-l X 0^0i^M0m/ 00 

— fe (fe — 1) Su^.'.'.'u^ x 0Mi"i ' ' • 0Mfc-2^-2 X 0" 0i^m0mA0a/0o 

+ x0m 1 , 1 ---0m^ (76) 

Factoring out 0q = —0a 0a we arrive at 

£fc (9^0)|«+l dimensions = "7 T~ ^ltV^fc+l X 0Ml I/ l ' ' ' 0Mfc-l l/ fc-l X 0Mfc + l0^fe+l 
with Q 0o=O (pA^A 



(77) 



n dimensions 



That is to say 



£ fc (590)| 

dimensions , , 

with a Q =O <PA0A 



n dimensions with 
varying a 0a 



(78) 



3.6 Legendre transformations 

The standard form for a Legendre transformation cf>, x < — > $, X is given by 



4>(x)+$(X) = '%2x a X a , (79) 



c*=l 



JMaO = ^Ua«0, z Q (X) = ^ffi> ee V Q $ . (80) 
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It follows that the Hessian matrices for and <E> are related by 

(ddcj))- 1 = (VV$) • (81) 
From this and the previous matrix identity (|37p it follows in n dimensions that 

' 1 £ k {ddcj>) = . j== . X — £ n _ k (W$) . (82) 



Vdct (0<90) k\ ^ ,/det (VV$) (n-fc)! 

That is to say, field equations for and $ are related by the Legendre tranform, and so are their solutions. 
The transformation gives a one-to-one local map between solutions of the nonlinear equations £ k {dd<j)) = 
and £ n -k (VV$) = 0, valid for all xorl such that the corresponding Hessian matrices are nonsingular, i.e. 
for all x or X such that det (000) ^ ^ det (VV$). 

This then is a general, implicit procedure for the construction of solutions to the equation £ k = given 
solutions to £ n -k = 0. In practice it is challenging to find tractable examples where the procedure can be 
fully realized. We will say more about solutions in §4. 

3.7 Legendre self-dual models 

The basic self-(anti)dual action consists of a pair of terms, 

A± = h Ak ± j^ky. An - k = J ^ (h. Zk - x {m) ± i^bfc)! £m W) ) • (83) 

Thus from (|43[) the first variation is 

SA± =-2js<f> £ k (500) ± ^ £ n - k (000) J ■ (84) 

This exhibits a classical self- (anti) duality for the resulting field equations, and their solutions, under the 
Legendre transformation (|79[) . Again in n dimensions, 

jmm 

- vmm (h £ * (VV *» ± <vv *>) ■ (85) 

In particular, for k — 1 this becomes 

' (£! (000) ± ^ (adj (000))) 



V'det (900) 

= (& (VV*) ± £i (adj (VV$))) , (86) 

Vdet (VV<J>) 



where we have made use of (|42|) . 

The consequences of this duality for quantized systems requires consideration of how the Legendre trans- 
formation directly affects the actions, A±, and not just the field equations. We consider this next. 

3.8 Legendre transformations of the action 

When M is taken to be the Hessian matrix, say, M = 000, then every £ k (000) is actually a double 
divergence, 

£ k (000) = 6$$*;::$* x Q1/3l Q2 ^ 2 • • • Qfc/3fc = d ai d Pl (a^X-ft,* x $ ' ' ' ^^*) • ( 87 ) 

Recall the Legendre transformation result 



gfeW)= (n- fc )!det(VV^) g "- fc ( VW) - (88) 
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Thus, using det = det (VV<1>), we also have 
J <j>, (x) 4>» (x) £ k (ddd>) d n x = J X.X, (n ! ! fc) , det( vv^) g "- fc (VV<I>) dGt (m) 



k\ 



(n 


-k)\ 




k\ 


(n 


- k)\ 




k\ 


(n 


— k)\ 




k\ 


(n 


-k)\ 



2 / KZ-2:: x * *°»p* ■ ■ ■ ^n-.p n -^x 



(89) 



where in the last step we used 

S;ZrZ = (n-m)S^Z Z . (90) 
The effect of the Legcndrc transformation is therefore 



Ak+i = J <t>n (x) ^ (x) £ k (ddcj)) d 
_ (fc + 1)! 



(n-k)l 
(fc + 1)! 



2 y$£„- fe -i(VV$) d n X (91) 



(n-k)\ 

After the usual integrations by parts, the latter Lagrangian has variation 

6 J C'-X-t x * *°»p> ■ ■ ■ *<*n- k Pn-j n x = {n-k)J (5$) x s^::.;::: $ Q2& • • • * an _ k0n _ k d n x , 

giving the expected equation of motion 

o = £„_*_! (w$) = ^i-Z~l ■ ■ ■ * an - kf > n - k ■ m 

It should be possible, therefore, to express the transformed Lagrangian in the standard form, upon integrating 
by parts: 

J s%:: v & x $ $„ lAl ■ ■ ■ $„ m x m d n x (93) 

= - J ^Ai-"-Am x $ ^a 2 ' ' ' ®v m \ m d n X (Ai integrated by parts) 

= - / OK:::*: - (m - 1) -t ) * (^^a.) <&„ 3 a 3 • • • ^ m x m d n X 

= -J ($A$a) (SZZZ >< ^A 2 ^ 3 A 3 ' ' ' $„ m A m ) d n X 

+ \ (m - 1) I x ($ Al ft^ ($ A * A )) $„ 3A3 • • • *„ mAra 

= - J ($A*A) (C-'.'.A™ X $„ 2 A 2 ^ 3 A 3 • • • ** m A J d n X 

- i (m - 1) /" ^AiAs-'.'.a" x (*a$a) $Ai^ 3 A 3 ' ' ' *^ m A m tP-X" (z/ 2 integrated by parts) 
So the re-expressed Lagrangian is simply 

J 6 v x \:Zl x $ • • • *^A m d"X = y $ £ m (VV$) d"X = -1 (m + 1) y (* A *A)f m -i (VV*) <TX . 

(94) 
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OK then, we have 
•Ak+i [4>] 

= J fa (x) fa (x) 8 k (ddcf)) d n x = j**^* 2 (~7^j (n — k — l + l) J ($a$a) £ n -k-i-i ( VV*) d n X 
= —, (fc t 1)! iM / (^$a)4-h-i(VV$) 

If we shift the index, this may be written more symmetrically. 

Thus we have established that the Legendre transform ([75)) gives directly a relation between the actions 
for the two theories: 

77 A [<t>] = 7^— ttt Ax-fc [*] Euclidean, (96) 
re! (n — k)\ 

provided boundary terms from integrating by parts may be discarded, where 

Ak [<i>] = j fa (x) fa (x) {dd<j>) d n x = ^ J <f> (x) £ k {ddcj>) d n x , (97) 

A n - k [*] = J $ M (X) (X) (VV$) d"X = - J k 2 + i /" $ (X) £„_ fc (VV#) . (98) 

But alas, the sign in (|9"6")l disagrees with the result from the explicit calculation of the n = 2, k = 1 case 
with Lorentz signature! Namely, 

'da^ d a (j) d 2 x = { V Q $ V Q $ d 2 X . (99) 



However, this discrepancy is due to the difference between the Euclidean and Lorentzian space identities 
for the product of two Levi-Civita symbols, namely, in n dimensions with Lorentz metric sign conventions 
(+,-,-,-,•••). 

e ai ... a „£ ft -^ = S^Zil Euclidean, (100) 
e ai ...a n e ft - A ' = (-if -1 S^Zil Lorentzian. (101) 

In Lorentz space then, the appropriate identity in n dimensions is 

1 f-D" 

TT A k [4>] = r TTT An-k * Lorentzian. (102) 

k\ (n — k)\ 
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3.9 Summary of results for Euclidean metrics 

Local relations: 



C 



k = (PaPa O k -l 



»ai^i0a 2 ^ 2 



*k-lPk-l 



£k = d p d a (S^l'.'.'.p k k X (f> a2l 3 2 ■ ■ ■ Qfc(3fc X 0^ B k = d a (^7fr--fc_ 2 2 X ^ift ' ' ' <?W_ 2 /3 fe _ 2 x P 0t0t) 

v k = d p d a {C£::Z-2 x ^ ■ ■ ■ ^ x ¥ 2 ) v k = ^\Z'-2:l&: x ■ ■ ■ x Mp> 

Note that £ k and T> k are double divergences. 



V k = § (1 - k) B k + i (1 + fc) £ fc V fe = P fe - 0£ fe 



Integrated relations: 



A k = j£k d n x 



Ak = T^fV k d n x 

JV k d n X : 

•Ak = jEpj j <t> &k d n x upon discarding both of the above boundary terms 



upon setting J B k d n x = 
J £ k d n x upon setting jT> k d n x = 



Constraint relation: 



Legendre relations: 



Trace relations: 



£k (990) |n+l dimensions 
with ch a d> a —0 



<t>\4>\ 



V fe+ j 



n dimensions with 
varying 4>\<p\ 



<t>(x)+$(x) = j: n a=1 x a x, 

d a 4> , x a (X) 



h £ k (990) 



d<S>(X) 
dX a 

1 



F- 4 * M ^ -JT^hy. A n-k [*] 



T ^n-fc (VV$) 



T fe = Tr (990) 



5 fe = Tr (9^90) (990)* 





/ 


71 


jfe-i 





• 





^ 








r 2 


Ti 


fc-2 • 


• 














r 3 


r 2 


Ti • 


• 










£/c = det 






















Tk-1 


7fe-3 


7fe-4 • 


• Ti 


2 











7fe-i 


7fe-2 


7fe-3 


• r 2 


Ti 


1 








r fe 


7fe-i 


7fe-2 


• r 3 


r 2 









/ 5 


k-i 





• 





^ 








Ti 


fc-2 • 


• 












s 2 


r 2 


Ti • 


• 










V k = det 


















Sk-3 


7fe-3 


7fe-4 


• Ti 


2 









S k -2 


T k -2 


7fe-3 


• r 2 


71 


1 






V s k -i 


7fe_i 


7fe- 2 • 


• r 3 


r 2 


Ti J 





Lorentz metric results depend on sign conventions, i.e. (+, —,—,■■■) or (— , +, +,•••), and are left to the 
reader to determine. 
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4 Classical solutions 



In this section we consider classical solutions of the field equations stemming from individual A^. We 
present and illustrate a variety of methods to find solutions of £k (dd(f>) = 0. 

As a warm-up, consider the first example beyond the standard free massless field, namely, £2 (ddcf>) = 0. 
For an extremely simple case, a "spherically symmetric" solution is valid almost everywhere in n dimensions: 

cj>{x) = { ^ a f a)l " 'f , (103) 

w \ \n(x a x a ) if n = 4 ' y J 

where a = 1, • • • , n is summed. However, in n-dimensional Minkowski space (except for n = 8, 12, • • • ) this 
solution obviously has an imaginary part outside the light-cone. So, in general, such functions with branch 
points give real solutions only on subspaces. We may think of the branch points as defining boundaries for 
the subspaces, with some particular boundary conditions. 
At the opposite extreme, another solution is 

4> (x) = ^xix 2 ■■■x n . (104) 

It is straightforward to check the equation of motion £2 — is satisfied. This solution has branch points for 
both Euclidean and Minkowski space. 

Next consider £ & = for general k. Although £k>i is nonlinear in 0, it is nevertheless still true that 
some plane waves are exact solutions. For "light-ray" plane waves, 

£ k [Aexp (iq a x a )] = (105) 

for constant A and q a , if q a q a — with A arbitrary. In this case, each of the terms in £k vanish separately. 
In fact, light-ray plane waves are only one among many possible solutions for which both <fi aa = and 
i'p'Pp — 0- That is to say, the general galileon equation £j~ = possesses a class of solutions given by the 
simultaneous solution of 

d 2 cj) n { dcj) 



, , - and — r - — r - )= . (106) 



The proof of this statement is elementary for £2 = 0, while for higher fc, the hierarchical construction in (1431) 
and the nature of the variation procedure guarantees that £k = will hold if both <j) a <j) a — and £k~i = 0, 
thereby establishing the general result^ 

In three dimensions it is only necessary to take the single constraint (j)a4>a — since a consequence of 
taking additional derivatives of this one constraint is <p aa = 0. This is not true in higher dimensions, but it 
does suggest another method for n = 3. 

4.1 Implicit solutions 

In three dimensions, choose four arbitrary functions f{u, v), g{u, v), h{u, v), k(u, v) constrained by the three 
relations 

xf(u, v) + yg{u, v) + zh(u, v) + k(u, v) — 4>(x, y, z), (107) 
xf(u, v) u + yg u (u, v) + zh u (u, v) + k u {u, v) = 0, (108) 
xf v (u, v) + yg v (u, v) + zh v (u, v) + k v (u, v) — 0. (109) 

Here subscripts denote partial differentiation with respect to u, v. Then the implicit solution of these 
equations for 4>(x, y, z) is a solution to the Monge- Ampere equation det {fypj) = in 3 dimensions. Here 

(j> x = f(u,v), (j> y = g{u,v), (j) z =g{u,v), (110) 

so the solution implies that there exists a functional relationship amongst {4> x , 4>y> 4>z)- This remark is 
enough to guarantee that the elimination of (it, v) will give a solution to Monge Ampere. To solve the 
galileon equation some further constraints are necessary. Suppose 

,<P y ) = Q(a,p) (111) 



4 It is interesting that this class of solutions is not given by the method of Legendre transformations, described above in §3.6 
and to be discussed further below. The Legendre method fails for this class because the second equation in d 106 1) implies the 
existence of a functional relation among the Legendre transformed variables so that det (dX a /dxfj) = 0. 
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defining a — <p x , /3 — <fi v . Then 

4>xz = Qa4>xx + Qfi4>xy , (112a) 
4>xz = Qa4>xy + Qp4>yy , (112b) 
4>zz = Qa4>xz + Qp4>yz , (112c) 
— (Qa) 4*xx 

+ 2Q a Q^ xy + (Q p ) 2 . (112d) 

Combining these relations we obtain 

{Myy ~ {<t>xy?) (l + Ql + Qp) , (113) 

so either the constraint is 

l + Ql + Q} = (114) 

or 

<t>xx4>yy - (fixy) 2 = . (115) 
In the case (|114[) this places a constraint upon the functional dependence of the derivatives. For example, 



if Q = i^J 4> x + 4>1 then the constraint is satisfied automatically. Then <j) x + (j)y + 4> 2 Z = 0, a constraint 

encountered earlier. In turn, this constrains the functions (f(u, v), g(u, v), h(u, v)) to satisfy f 2 +g 2 +h 2 = 0. 
Another possibility is Q = ism(cf> x ) + icos(4> y ). 



In the other case where (| 1 15|) holds it is not sufficient to impose just this to solve the galileon equations, 
but we also require that (j) z = Q{4> x , 4> y ), or else demand that all leading subdcterminants vanish, i.e. 

4>xx4>yy - (4>xy) 2 = 0, 4>yy4>zz ~ {4>yz) 2 = 0, 4> Z z4>xx ~ {4>zx) 2 = . (116) 
A solution of this type is given by <p(x, y, z) = ^Jxyz. 

4.2 Envelope method 

The envelope method also gives solutions. To be explicit, in 3D, the method is to take 

xf (u) +yg(u) + zh(u) = <i>(x,y,z) , (117) 
xf'(u) + yg'(u) + zh'(u)=0, (118) 

and then choose various functions /, g, and h to determine u (x, y, z). For example, inserting / (u) = u, 
g (u) = 5 it 2 and h (u) = -|u 3 into the second of these two equations gives x + yu + zu 2 = 0, whose solutions 
are 

u (x, Vi z ) = y z (-V ± vV - 4xz) . (119) 
Therefore the corresponding solutions for <f> are 

4> (x, y , z) = -^22 {~y 3 + 6x y z ± (y 2 ~ 4xz ) ^y 2 ~ 4xz ) ■ ( 12 °) 

Again, these functions have a branch point, as well as a pole, so they are not real solutions unless y 2 > 4xz 
and z ^ 0. For these examples, it is again straightforward to check the equation of motion is satisfied. 

This procedure can be extended if all three derivatives are functionally related, and the auxiliary functions 
depend only upon one function: 

xf(u) + yg{u) + zh(u) + u)(u) = <f>(x, y, z) , (121) 
xf u (u) + yg u {u) + zh u {u) + uj u (u) = 0. (122) 

The first of these equations possesses a tantalizing similarity to the Legendre transform. But here the 
procedure is to solve the second equation for u in terms of x, y, z which results in a solution of £2 = when 
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F(u) 


= q(u)f(u) 


G(u) 


= q{u)g(u) 


H(u) 


= q(u)h(u) 




= q{u)ui(u) 



inserted into the first equation. However, having obtained this solution for u, a huge class of solutions may 
be obtained by replacing f{u), g(u), h(u), and lo(u) by 

(u)f(u)du, (123) 

(u)g(u)du , (124) 

{u)h(u)du , (125) 

{u)uj{u)du . (126) 

The equation to determine u is still the same as before, so a new solution is generated. This may be checked 
on the specific example, taking q{u) = u. 

Another example 

Take f(u) = u, g(u) = u 2 , h(u) — u 3 . Solving the second of these equations for u, 



u = l/3 - y+Vy2 ^ , (127) 

z 

and 

( — V + \fy 2 — 3zx\ (6zx— y 2 + y\J y 2 — 3 zx) 

0=1/27-^ J -± — 5 '- . (128) 

z l 

This is a solution to the n — 3 case. Note this is of weight one, so it is also a solution of the Monge- Ampere 
equation. 

4.3 Power law solutions in other dimensions and Legendre equivalences 

In general, the equation £ k (dd(f>) = is homogeneous in 0, of degree k, and therefore the overall normalization 
of any solution is not determined. Consider again spherically symmetric solutions as given by a power ansatz: 

<t>(x) = {x a x a f . (129) 

In n dimensions, for this ansatz, the products and traces of ddifi are: 

(dd4>) k ^ = {2p) k {x a x a 5^ u + ((2p- l) fc - l) Xf^xS) {x a x a ) kp ~ k ~ 1 , 

(00$J M = {2pf (n-l + (2p- if) (x a x a ) kp - k . (130) 

Inserting these traces into (|2l)f and evaluating (|2"5j) we find, for example, 

£x = 2p (n - 2 + 2p) (i a i tt ) rl = p = 1 - n/2 , (131) 
£ 2 = Ap 2 (n - 1) (n + 4p - 4) (x a x a ) 2p ~ 2 = p = 1 - n/4 , (132) 
£ 3 = 8p 3 (n - 1) (n - 2) (n + 6p- 6) (a;^) 3 ^ 3 = p = 1 - n/6 , (133) 

Thus the power p required for a solution is determined, as indicated. 
For other levels in the hierarchy, 



£k W)l, (a0 =(»«,»«,)' = (x a x a ) kp - k (n + 2kp - 2k) (2p) k ^ 2) f ' 
The condition on p for the ansatz (|129p to be a solution of £ k = is therefore 



(134) 



By taking a limit, a nontrivial solution for n = 2k is easily found to be In (x a x a ). 
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So, in n dimensions the fcth equation of motion of the hierarchy is solved by 

(136) 

In (x a x a ) if n = 2k 

Moreover, under the Legendre tranformation, the ansatz solution for level k is mapped into the ansatz 
solution for level n — k. 

In particular, when n = k + 1 the ansatz is mapped into an harmonic function by the Legendre transfor- 
mation!^ In n dimensions the spherically symmetric harmonic solution is given by = V^V^ (X a X a ) 2 , 
so this last statement is equivalent to 

(i a x tt )^ = (x a x a )^ ~ (X a X a )^ = (X a X a )^ (139) 

Legendre for n—k+1 

Note the effect of the transformation is just to replace k i — > i in the exponent, thereby changing the scaling 
properties of the solution: 

(XpVp-ljiXaXa) 2 -^ 1 ={l-n){X a X a ) i ^ L , (140) 
(xpdp - 1) (x a x a )^^ = (t^) t^ 1 ") 5 ^ ■ ( 141 ) 
Let's go through the details for the n = k + 1 case. Under the Legendre transformation: 

xpXp = 4> (x) + $ (X) , (142a) 

3d) (x) fc— i — fe— i 

Xp = = — (x a x a ) 2 " = ^-xp (x a x a ) 2 <= , (142b) 

OXp OXf) K 

x p Xp = h^L (x a x a )^ =^4>{x) , (142c) 

Z[X) = -\ <P(x) , (142d) 

X fj Xp = (^i) 2 (x a x a )^ , x a x a = (^) 2fe {XpX y k , (142c) 

(x a x a )^ = (^i)^ 1 {XpXp) 1 ^ . (142f) 

So then, an harmonic function of X is indeed the result of transforming (x a x a ) 2k , for dimension n = k + 1. 
Including a convenient normalization, 

ct>{x)=k k {x a x a )^ ^{X) = -(k-l) k - 1 (Xi 3 X p ) 1 -^ . (143) 

Legendre for n—k+1 

Note that this procedure could be reversed, starting from the harmonic solution. Thus there is a local 
one-to-one map between harmonic functions $ (X) and solutions of the nonlinear equation E n -i (dd(f>) = 0, 
in n dimensions, valid so long as < |det (VV$)| < oo. 

4.4 Self-Dual Solutions 

Another approach to the solution provides a class of self dual solutions; i.e solutions both to the original 
equation, and to the Legendre transformed equation (in the same variables). Suppose we impose the ansatz 

xM4>) = 1. (144) 



5 Two such examples for solutions of (dd(f>) = and Si (VV<J>) = in three dimensions are 

<!> = XYZ , <j> = ^xyz , (137) 

and 

* = 7WTWTW> +=<*W + *) 1 '*. das) 

For another, take <t> = Z (Z 2 + Y 2 — 4X 2 ), etc. We leave it as an exercise for the reader to find the corresponding (f>. 
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(Here the subscript fi is a vector index, not a derivative.) Then 

(145) 
(146) 



d± = __K L _ 

(VM + Kvp , TOg) 



Inserting this expression into £2 we obtain 



(E 



(EWE W) - E W) ■ ( 147 ) 



By construction solutions of this type will inevitably be complex in Euclidean space, but may be real in 
Minkowski space. This is clearly zero if the constraint E VfiVp, = is imposed. Moreover, it is easy to see 
that this solution also solves V 2 = 0, the Legendre transform of £2 in 3 dimensions. Another remarkable 
property jTB] of this solution is that if you replace <f> by any function of (f>, it remains a solution of these 
equations! Indeed, this class of solutions is universal, as it solves all equations £k = in the appropriate 
dimension, as every term contains at least one factor of the form E V^V^ or E V/jVu) which both vanish by 
the constraint. Furthermore this extends to covariant equations which also include first, as well as second 
derivatives; thus the Lagrangians themselves vanish on this class of solutions. 

Example 

As a simple example consider the equation 

x^ + iy^Jl + cj) 2 + z = 1 (148) 
whose coefficients satisfy the constraint. Solving for 4> we obtain 



(149) 



±{x — zx) + -J — y 4 — y 2 z 2 + 2 y 2 z — y 2 — x 2 y 2 

( P= — 2 ■ 

y + x z 

These solutions may be verified to satisfy £\ = 0, £2 = 0, and £3 = 0. 

5 Mixtures 

In this section we consider solutions of the field equations for systems governed by linear combinations of the 
Ak for different k. 

5.1 Symmetric spacetime solutions 

Consider the modified field equation that follows from Ai — kA2- 

<t>\\ = k (<p vv <p\\ — 4>Xv<P\v) ■ (15°) 
Make a spherically symmetric spacetime ansatz: 

<f> = f(a) , a = x x x\ . (151) 

Then 

<j>xx = d x (2a*/') = 2nf + 4x x xxf" , fax = d v (2x x f) = 28 vX f + 4x v xxf" , (152) 

4>xxfa, - faxfax = (2nf + Axxxxf"? ~ (25 vX f + 4x v x x f") {28 vX f + 4x v x x f") 

= An (n - 1) (f'f + 16 (n - 1) x x xxf'f" , (153) 

faxfa* - faxfax = 4 (n - 1) (/') (nf + 4af") , (154) 
and the field equation (|150[) becomes 

2nf + 4a f" = 4n{n- 1) (/') (nf + 4a f") . (155) 
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This is again a first order differential equation for g = f: 

2 , 2«(n-l) 2 , ir „, 

-erg' = -g+ v - ' - n 2 . (156 

n (4k (n- l)g - l) y v ' 

Therefore 

in K 5+ Rib)))=-? in ^ ^ 

and we have 

9 ( 9 + ji^) ) ^ ^- /2 where C = 5l + ^i— y ) , . 9l .,(. = !). (158) 



That is to say, 



1 ( 1-a/i + 16k 2 (n-l) 2 Ca-«/ 2 if 4«n-l ?1 <l 

.9( CT ) = -r-^ TTS / 5 • ( 159 ) 

i n ~ 1 ) { 1 + y'l + 16k 2 (n - l) 2 Ca-"/2 if 4k (n - 1) 9l > 1 



4k 

One more integration gives / = J g: 



i» = J>o) + I g(p)dp. (160) 

•J an 



'on 

For convenience, let us take <to = 1. The result of the integral is then 

l + m K 2 (n-l) 2 Ca- n / 2 d<j = <j 2 F 1 (--, --:l--;-16K 2 (n-l) 2 Ca-^A , (161) 

\ 2 n n ) 

where iF\ is the usual Gauss hypergeometric function, 

qfrz a (l + q)6(l + b)z 2 a(l + a )(2 + a)b(l + b)(2 + b)z 3 
2 F 1 (a,6,c,,) = l + — + c(1 + c)2 , + C( i + C )( 2 + C )3! + " " (162) 

The final result for the spherically symmetric spacetime solution (recall a = x\x\) is then 

, ±i f 2 ^(-|,-f;i-f;-i6K 2 (n-i) 2 c) \ 

4« (n - 1) ^ _ ff aFl (_ i , - 2 ; 1 - 2; -16k 2 (n - I) 2 Ca~^) ) ' 1 J 

where the ±1 choice is made depending on whether 4k (n — 1) gi ^ 1, thereby giving various values for 

— 16k 2 (n — 1) C — (2 — 4k (n — l)<?i) x 4k (n — At the critical value 4k (n — 1) #1 = 1 we have 

— 16k 2 (n — l) 2 C = 1. (Note: The series expansion is not valid for the hypergeometric function if 
16k 2 (n - l) 2 Ca~^ n > 1.) With z = 4k (n — 1) <?i we have —16k 2 (n — l) 2 C = z (2 — z), and we note 
for z (2 — z) = — 1, the solutions are: z — 1 ± \/2. 

5.2 Static, spherically symmetric solutions 

For example, for the free field in 4 spacetime dimensions, S\ = 0, the static, spherically symmetric solutions 
of V 2 (p (r) — in 3 space dimensions are of course 

<£i (r) = C + — . (164) 
r 

For the next step up in the hierarchy, £2 = 0, the static, spherically symmetric solutions in 3 space dimensions 
are 

02 (r) = C + C 1 ^ . (165) 

These two solutions are Legcndrc duals in 3D space (but not in 1+3 spacetime). So, what happens if we 
mix them up? 
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For example, take 

£1 = \£-2 ■ 

The static, spherically symmetric solutions in this case satisfy 



^2 d r (r 2 d r 4>) = ^d r (r (d r (j))' 



which has an immediate first integral, and solution, 

d = r 2 d r <\> + 2Ar {d r 



Integrating this gives 



where the length scale is related to the previous first integral by 

R = (8ACi) 1/3 . 



117 r 3 

2' 6' 6'~B3 



(166) 
(167) 

(168) 
(169) 

(170) 
(171) 



For small r the solution (|170j) behaves like 02 (r) in (| 165[) . and therefore it is not singular at the origin. 
On the other hand, for large r, upon taking the upper + sign in (|170|) . the solution behaves like 4>\ (r) in 
(|164l) . while the lower — sign choice in (|170[) gives a solution that grows like r 2 for large r. 

Taking the upper sign in (11701) . 



R 2 



\ V 2 



where 



3r(i/3)r(7/6) 



^K:^ + itx{ 4 \Ir-v +0 ((r) 

Hr) m + f x 3r(i/3)r ( 7/ 6 ) _g + //g 



4. 2065 • • • . Taking the upper sign, the graph of 8 A (<f> (r) ~ (f)(0)) /R 2 follows. 



(172) 
(173) 




f£ ((/) (r) - (f> (0)) versus r/R, for R > 0. 



Another branch, obtained by taking the lower sign in (|170j) . is not so well-behaved for large 
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Two branches of the mixed solution, J=£ (0 (r) — <f> (0)), versus r/i?, for R > 0. 
5.3 Energy considerations 

A symmetric energy-momentum tensor for the mixed £ \ — £2 model is 

d^O^ = E m <f>, , (175) 

£ W\ — §aa - A (0QQ0/3/3 - (fra^ap) , (176) 

where the equation of motion is £ [4>] — 0. For static, spherically symmetric </>, the energy density is 

Ooo = \ {dA (r)f - ^d r {d r <f> (r)) 3 , (177) 

and the total energy is 



E = Att 6oo r 2 dr . (178) 
Jo 

This is finite for the bounded static solution that goes like 1/r for large r, but it is not finite for the solution 
that goes like r 2 . For the finite case, 



^=^(-i+V 1+ ?j ' (179) 

where R 3 = 8XCi gives a length scale set by the first integral of the static equation. After some playing 
around, we find 



nR 5 [°° , / / l\ 1 2^tt 3 R 5 R 5 , N 

This is true for R > 0, but actually, it is also of interest to consider cases where R < 0. 
The bounded solution for R < is real for r > \R\, as obtained by integrating 



23 



Two branches of the mixed solution, (0 (r) — cj) {\R\)) 1 versus rj \R\, for R < 0. 

The energy density has the same form as before, but now the total energy outside the singularity at 
r = \R\ is 



E, 



>\R\ 



6A 2 



-1 + 1 1 



1 



ds + r 2 (9 r 0) 3 



r=|Ji| 



237T 3 3 

7T 



180 V(r(f)) 3 4 y A 2 



/?l = 0.097 037 ^ 



(182) 



5.4 Perturbative scattering 

Consider p + q — > p' + q' for the L\ + XC2 model, perturbatively on-shell, i.e. p 2 = q 2 = p' 2 = q' 2 = 0. The 
lowest-order scattering amplitude is 



M = \x 2 ( S 3 + t 3 + U 3 ) \ u= _ s _ f = \\ 2 stu\ u _ s _ t = -±#st (s + t) 



(183) 



In the CM frame, in terms of the incident energy and scattering angle, p = (E, and p ■ p' = E 2 cos 9, 
we have s = AE 2 , t = —2E 2 (1 — cos#), and so 

M CM = 12A 2 £ 6 sin 2 6 = 12A 2 £ 6 (l - cos 2 6) = 8X 2 E 6 (P (cos 9) - P 2 (cos 9)) . 

By the usual rules for the differential cross section in 4D, we then have 



da 

dh 



1 1 IMcmI 9A 4 £ 10 sin 4 (9 



CM 



64tt 2 2! AE 2 



32tt 2 



and total cross section 



4 P 10 



I CM 



3A 4 £ 
5tt 



(184) 



(185) 



Note that [A] — 1/m 3 in 4D. This approximation for a obviously exceeds the Froissart bound (cx In 2 E) as 
the energy increases. 
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6 Effects of cf) Q[4>] self-couplings 



In this section, we consider galileon theories with an additional self-coupling of the fields to the trace of their 
own energy-momentum tensor. We explore the classical features of one such model, in flat 4D spacetime, 
with emphasis on solutions that are scalar analogues of gravitational geons. We discuss the stability of these 
scalar geons, and some of their possible signatures, including shock fronts. 

For the simplest example, the galileon field is usually coupled to all other matter through the trace of 
the energy-momentum tensor, Q( mattG1 ). But surely, in a self-consistent theory the galileon should also be 
coupled to its own energy-momentum trace, even in the flat spacetime limit. Some consequences of this 
additional self-coupling are considered in this section, based on work published in [8] . 

Recall the action for the lowest non-trivial member of the galileon hierarchy, 

A * = \ j taMpp d n x , (186) 

where (j) is the scalar galileon field, <p a = d(f>{x) /dx a , etc., and where repeated indices are summed using 
the Lorentz metric <5^„ = diag (1, — 1, — 1, • • • ). 

6.1 A non- vanishing trace 

As discussed above, including in A<i a minimal coupling to a background spacetime metric yields a symmetric 
energy- momentum tensor, which becomes in the flat-space limit: 

®JS = <t>H<P» ( t>aa ~ <$>a.§av<l>n ~ 4>a4> at i4>v + & nAa^pfyafi ■ (187) 

This is seen to be conserved, 

0„e$ - <f>„ £2 [<t>] , (188) 

upon using the field equation that follows from locally extremizing A2, = bAijb'fy = — £2 [<t>], where 

£2 [4>] = <f>aa<t>/30 - 4>ap4>al3 ■ (189) 

(2) 

But, as previously noted, this Q^J is not traceless. Consequently, the usual form of the scale current, 

( 2) ______ 

x a Q ai i, is not conserved [26] . On the other hand, the action (|186p is homogeneous in <f> and its derivatives, 
and is clearly invariant under the scale transformations x — > sx and 4>(x) — > s <i4 ~ n ^ 3 (j)(sx). Hence the 
corresponding Noether current must be conserved. This current is easily found, especially for n — 4, so let 
us restrict our attention to four spacetime dimensions in the following. 
In that case the trace is obviously a total divergence: 

e< 2 > = ve{fi - d a faa^) . (190) 

That is to say, for n = 4 the virial is the trilinear V a = 4> a 4>i34>i3- So a conserved scale current is given by 
the combination, 

s» = s a ew - . (191) 

Interestingly, this virial is not a divergence modulo a conserved current, so this model is not conformally 
invariant despite being scale invariant. Be that as it may, it is not our principal concern here. 

Our interest here is that the nonzero trace suggests an additional interaction where <f> couples directly 
to its own O^ 2 ). This is similar to coupling a conventional massive scalar to the trace of its own energy- 
momentum tensor [21 . In that previously considered example, however, the consistent coupling of the field 
to its trace required an iteration to all orders in the coupling. Upon summing the iteration and making 
a field redefinition, the Nambu-Goldstone model emerged. But, for the simplest galileon model in four 
spacetime dimensions, (|186|) . a consistent coupling of field and trace is much easier to implement. No 
iteration is required. The first-order coupling alone is consistent, after integrating by parts and ignoring 
boundary contributions, so 

-J / (t> d a {4> a (f>p(f>fs) d^X — j / (t> a 4>a4>l3(t>li ^X . (192) 



3 Also note that A2 follows from coupling <j> to the trace of the manifestly chargeless tensor (d^du — <5 M i/9 a c9a) <f>s<t 
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(Similar quadrilinear terms have appeared previously in [TU1 [H] , only multiplied there by scalar curvature R 
so that they would drop out in the flat spacetime limit that we consider.) Consistency follows because (|192[) 
gives an additional contribution to the energy-momentum tensor which is traceless, in 4D spacetime: 



6$ = Mu<f> a 4> a - \S^<t>*<t>cc<t>p<t>f> , © (3) = • (193) 

Of course, coupling <fi to its own trace may impact the Vainstein mechanism |36) by changing the effective 
coupling of 0( mattcr ) to both backgrounds and fluctuations in (f>. We leave this as an exercise for the reader. 

6.2 A model with additional quartic self-coupling 

Based on these elementary observations, we consider a model with action 




(194) 



where for the Lagrangian L we take a mixture of three terms: the standard bilinear, the trilinear galileon, 
and its corresponding quadrilinear trace-coupling. The quadrilinear is reminiscent of the Skyrme term in 
nonlinear a models 33J although here the topology would appear to be always trivial. 

The second and third terms in A are logically connected, as we have indicated. But why include in A 
the standard bilinear term? The reasons for including this term are to soften the behavior of solutions at 
large distances, as will be evident below, and also to satisfy Derrick's criterion for classical stability under 
the rescaling of x. Without the bilinear term in L the energy within a spatial volume would be neutrally 
stable under a uniform rescaling of x, and therefore able to disperse [TT1 H4] . 

Similarly, for positive k, the last term in A ensures the energy density of static solutions is always bounded 
below under a rescaling of the field 0, a feature that would not be true if k = but A ^ 0. So, we only 
consider k > in the following. But before discussing the complete 6^ for the model, we note that we did 
not include in A a term coupling (f> to the trace of the energy-momentum due to the standard bilinear term, 
namely, J (f>Q^d 4 x : where 

= Kh> ~ 3<W«&« . 0(1) = -Ma ■ (195) 

We have omitted such an additional term in A solely as a matter of taste, thereby ensuring that L is invariant 
under constant shifts of the held. Among other things, this greatly simplifies the task of finding solutions 
to the equations of motion. 

The field equation of motion for the model is = 6 Aj '5(f> = —£ [<f>], where 

£ [(f)] EE (f> aa - A ((f> aa <f>/3f3 - (f>al3<f>ap) ~ K ((f>a4> '/3<P 73 ) a ■ (196) 

As expected, this field equation is second-order, albeit nonlinear. Also note, under a rescaling of both x 
and <f>, nonzero parameters A and k can be scaled out of the equation. Define 




(197) 



Then the field equation for ip (z) becomes 

Ipaa ~ (4>aalf>l3l3 ~ ^a^ap) ~ (^ ) a'0/3#) Q = , (198) 

where ip a = dip (z) /dz a , etc. In effect then, if both A and k do not vanish, it is only necessary to solve the 
model's field equation for A = n = 1. 

6.3 Static solutions 

For static, spherically symmetric solutions, (f> = (f>(r), the field equation of motion becomes 




(199) 



where (f>' = d(f>/dr. This is immediately integrated once to obtain a cubic equation, 

rV + 2Ar {4>'f + nr 2 {(f>'f = C , (200) 
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where C is the constant of integration. Now, without loss of generality (cf. (|197[) and Q198P ) we may choose 
A > 0. Then, if C = 0, either <f>' vanishes, or else there are two solutions that are real only within a finite 
sphere of radius r = y/ A 2 /k. These two "interior" solutions are given exactly by 

<f>' ± = -— (x ± VA 2 - r 2 n) . (201) 



Note that these solutions always have <fi' < within the finite sphere. 

Otherwise, if C ^ 0, then examination of the cubic equation for small and large \4>'\ determines the 
asymptotic behavior of 4>' for large and small r. In particular, there is only one type of asymptotic behavior 
for large r: 

C 

4>' ~ — for either sign of C . (202) 

r— >oo T 

However, there are two types of behavior for large \4>'\, corresponding to small r. Either 

2Xr i + o(±)) (203) 



4>'k 

provided <p' < 0, but with either sign of C; or else 

r = Mfi + (^)) (204) 

provided C > 0, but with either sign of <fr' . The corresponding real solutions behave as 

— 2A 

for either sign of C, or (205) 



r->0 Kr 



o ± V 2~^ provided C > ' ( 206 ) 

Comparison of the small r behavior to the large r asymptotics shows that in half these cases the solutions 
would require zeroes to be real and continuous for all r. But such zeroes do not occur. Instead, half of the 
cases provide real solutions only over a finite interval of r, somewhat similar to the C = solutions in (1201[) . 
but not so easily expressed, analytically. 

The solutions which are real for all r > boil down to two cases, with small and large r behavior given 
by either 



C , „ C 
2Ar 

or else 



and 4> _ ~_ for C > 0, (207) 



— 2A C 
0' an d ft ~ — forC< 0. (208) 

i — >0 Kr r— >oo r 

From further inspection of the cubic equation to determine the behavior of 0' for intermediate values of r, 
when C > it turns out that cf>' is a single- valued, positive function for all r > 0, joining smoothly with 
the asymptotic behaviors given in (|207l) . However, it also turns out there is an additional complication 
when C < 0. In this case there is a critical value (k 3 / 2 /A 2 ) C cr iti ca i = — 4\/3/27 ~ —0.2566 such that, if 
C < C cr iticai then <fi' is a single- valued, negative function for all r > 0, while if C cr iticai < C < then <j>' 
is triple- valued for an open interval in r > 0. It is not completely clear to us what physics underlies this 
multivalued-ness for some negative C . But in any case, when C < it is also true that <fi' joins smoothly 
with the asymptotic behaviors given in (|208|) . All this is illustrated in Figures 1 and 2, for A = K = 1. 

A test particle coupled by 00( mattcr ) to any of these galileon field configurations would see an effective 
potential which is not 1/r, for intermediate and small r. Therefore its orbit would show deviations from 
the usual Kepler laws, including precession that is possibly at variance with the predictions of conventional 
general relativity. It would be interesting to search for such effects, say, by considering stars orbiting around 
the galactic center [28l [27] . 
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Figure 1: ip' (r) for C = -l/2 N , with TV = 6, 5, 4, 3, 2, 1, from left to right, respectively. The thin black 
curve is a union of the two C = solutions in (|201j) . 




■0' (r) for C = +1/4 N , with iV = 0, 1, 2, 3 for top to bottom curves, respectively. 

For the solutions described by (|207[) and (|208[) . the total energy outside any large radius is obviously 
finite for both C > and C < 0. And if C > 0, the total energy within a small sphere surrounding the 
origin is also manifestly finite. But if C < the energy within that same small sphere could be infinite 
unless there is a cancellation between the galileon term and the trace interaction term. Remarkably, this 
cancellation does occurQ So both C > and C < static solutions for the model have finite total energy. 

6.4 Energy considerations again 

Complete information about the distribution of energy is provided by the model's energy-momentum tensor, 

%u = e$ - Ae$ - «e$ . (209) 

7 For C < 0, to see cancellation between the individually divergent galileon and trace interaction energies for small r requires 
leading and next-to-leading terms in the expansion: <f>' ~ + O ( r )- 
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As expected, this is conserved, given the field equation £[</>] = 0, since 



a M e^ = 4> u £ [</>] . (210) 

The energy density for static solutions differs from the canonical energy density for such solutions (namely, 
—L) by a total spatial divergence that arises from the galileon term: 

©oo = - iUtic - • ((V<P) 2 ^A . (211) 



This divergence will not contribute to the total energy for fields such that linv-^oo ((/>/ lnr) exists. Assuming 
that is the case, Derrick's scaling argument for static, finite energy solutions of the equations of motion |llj 
shows the energy is just twice that due to the bilinear Q^j . Thus, 



E = I Onn d 3 r 



J (^</>) 2 d 3 r . (212) 

For the spherically symmetric static solutions of (|200|) . this becomes an expression of the energy as a 
function of the parameters and the constant of integration C : 

/>oo 

E [A, k, C] = 4tt / (<f>'f r 2 dr . (213) 
Jo 

Again without loss of generality, consider A = K = 1. Then for either C > or for C < C cr iti C ai < 0, change 
integration variables from r to s = cf>' to findU 

E(C%0)=I(\C\)t(\C\ + ±v) , (214) 

I (C > 0) s \ r ds , (215) 

2 Jo (s 2 + l) 4 ^ s 4 + s(s 2 + 1)c 

where the numerator of the integrand is an eighth-order polynomial in s, namely, 

P (s, C) = 8s 8 + 12C*s 7 + (3C* 2 - 8) s 6 + 8C*s 5 + 7C 2 s 4 - 4Cs 3 + 5C 2 s 2 + C 2 . (216) 

Thus, / (C) is an elliptic integral. But rather than express the final result in terms of standard functions, 
it suffices here just to plot E (C), in Figure 3. Note that E increases monotonically with \C\. 

For other values of A and k with the constant of integration C specified as in (|200|) . the energy of the 
solution is given in terms of the function defined by (|214I215|) : 

E[\,k,C\ = (a 3 /« 5/2 ) e[k 3 ' 2 C/\ 2 ) . (217) 

The energy curves indicate double degeneracy in E, for different values of |C|, when E [A, k, C] > 7tA 3 /k 5 / 2 . 
Also, for a given |C| the negative C solutions are higher in energy, with E [A, k, — |C|] — E[X,k, |C|] = 

7tA 3 /k 5 / 2 +2 |C| X/k. Or at least this is true for all \C\ > |C cri ticai| in which case E [A, k, C] > -^E ^^^-Ccriticai 

3.7396 A 3 /k 5 / 2 [?]. 



6.5 Scalar geons and a shock-front conjecture 

Finite energy classical solutions of gravity-like theories bring to mind the "geons" proposed long ago by 
Wheeler 39,. These were envisioned in their purest form as distributions of only gravitational energy 
held together solely by gravitational interaction. Combinations of electromagnetic energy and gravity were 
also considered, as were systems containing neutrinos. Wheeler argued that such configurations would 
be relatively stable, if they existed, but would eventually dissipate due to a variety of both classical and 
quantum effects, including light-light scattering, as well as production and absorption of quanta. While 
plausible distributions were sketched, and decay rates were estimated, exact classical solutions were not 
found. 

8 The multivalued behavior of any solution for C cr ; tical < C < makes the determination of the total energy ambiguous, at 
best, for these cases. This is an unresolved issue. 
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C 
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Figure 2: E (±C) versus C > as lower/upper curves (the horizontal line is E (C cr iti ca i) ~ 3.7396). 



The same mechanisms would seem to apply to any hypothetical classical galileon distributions such as 
those discussed here, the main difference being that analytic spherically symmetric solutions might still be 
obtainable even if conventional gravitational effects were included. Perhaps these gravitational effects would 
not alter the qualitative features of the static pure <j> configurations given above. Should they really exist, 
presumably these galileon geons could also be dissipated by various classical and quantum effects. All this 
is far beyond our current abilities and the scope of this paper, of course, but the general ideas suggest some 
interesting possibilities. 

Whatever the cause, if the configuration's energy loss were gradual, as a first step it might suffice to model 
the time-dependent system quasi-statically, as a continuous flow from one static solution to another. That 
is to say, perhaps a good approximation would be to take C (t), with \C\ and E (C) decreasing monotonically 
with time. For the positive C case, this would be more or less uneventful as the whole configuration would 
just slowly disappear without any abrupt changes. But for the negative C case, as t increased Critical would 
be reached, beyond which the solution would begin to fold over, exhibiting the multivalued features shown 
in the Figure. But this is just the usual picture for the formation of a shock front. These particular galileon 
shocks would implode, converging towards the origin, as shown Iherel We believe this is a plausible scenario 
and a reasonable physical interpretation of the model's multivalued solutions. Moreover, it would seem to 
provide a signature for their existence. 

As is clear from the Figure, the shock front would form when dip' /dr — oo. For the C < static solutions 
of (|200j) it is not difficult to determine the locus of such singular points. It is given by the intersection of the 
solutions, for various C, and the curve (l -I- 3k</>' 2 ) r — AX(f>' . As usual for singular points in the development 
of a shock, almost certainly there is some physics missing from the equations. Since 4>" is large, the obvious 
modification would be to include higher derivative terms in the action, which is tantamount to attempting 
an ultraviolet completion of the model. This is an open question. Perhaps higher terms in the galileon 
hierarchy would be natural candidates to be included. 

6.6 Comparison to the self-dual model 

To get a handle on such terms, and for purposes of comparison to the model in (|194j) . consider briefly 
another model somewhat similar in form, but whose Lagrangian consists only of terms taken from the galileon 
hierarchy, without any coupling to 0. After rescaling the field and coordinates to achieve a standard form, 
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this alternate model may be defined by 

^self-dual [tp] = J {\lll a Tpa ~ \lpai>ai>l3l3 

+ J2tpa4>a (V , /3/3'077 ~ IpPf^Pf)) ^ X ' ( 218 ) 

The difference with (| 194|) lies in the last term, which is quadrilinear in the field, as before, but now has two 
fields with second derivatives. 

As the name suggests, this model is self-dual, in the following sense: The action retains its form under 
a Legendre transformation [TH] (also see to a new field \& and new coordinates X, as defined by: 

V> (x) + * (X) = x a X a . (219) 

Thus j4scif-duai [ip] — ^4scif-duai provided integrations by parts give no surface contributions. This identity 
suggests that there are interesting properties for the quantized model, such as its ultraviolet behavior, but 
that is outside the scope of the present discussion. 

Here it suffices to compare the classical physics following from (|218|) with that following from ()194|) . 
Upon integrating once the classical equations of motion for static, spherically symmetric solutions of the 
field equations for (|218l) . the result is again a cubic equation, 

rV+r(^) a + jW 3 = C, (220) 

but the (ip 1 ) 3 term is no longer weighted by r 2 as it was in (|200[) . Thus the small and large r behaviors are 
now given by 

tp' ~ (3C) 1/3 and i>' ~ (221) 

) — >0 r— >oo T 

for either sign of the constant of integration, C . These static solutions have finite total energy for either 
sign of C, as before, only now -0' is always bounded. Moreover, upon inspection of the behavior of ip' 
for intermediate r, and various C, unlike the previous model the solutions are now always single- valued for 
either C > or C < 0. Thus there are no multivalued solutions like those shown in the previous figure for 
various C < 0. However, each of the C < static solutions now has a single point for which dip' /dr = oo, 

1/3 

namely, r = (3 |C|) ' . So there is still a reason to expect the existence of shock fronts for quasi-static 
time-dependent fields in this alternate model. Finally, again for C < 0, to have (/>' real for all r > 0, it 

1/3 

is necessary to join together "interior" and "exterior" solutions at r = (3 |C| /2) ' . These features are 
illustrated in the following figure. 
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Static solutions of the self-dual model for C > upper half-plane, and for C < 0, lower half plane. The 
solutions r± — y 1 ± ffy* are snown m orange/green. 

In these graphs, the solutions of (|220l) are shown for both physical r > and unphysical r < to display 
some symmetry relations between the C > cases and the interior and exterior solutions for C < 0. The 
straight lines, in gray, are the loci of points where the solutions have zero and infinite slopes, for different 
values of C. 

It remains to investigate the stability of these spherically symmetric solutions under perturbations, es- 
pecially to check for the existence of superluminal modes, along the lines of |23j . Evidently, superluminal 
modes are a possible feature for models of this type. 
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7 General relativistic effects 



In this section, the simple trace- coupled Galileon model of the previous section is coupled minimally to 
gravitation ( GR ) and shown to admit spherically symmetric static solutions with naked spacetime curvature 
singularities. 

In the previous section, based on [5] , the effects of coupling a Galileon to its own energy- momentum trace 
were considered in the flat spacetime limit. Here, general relativistic effects are taken into consideration and 
additional features of this same model are explored in curved spacetime 110) . Such features have been 
explored in the literature (see [7J, and for a related class of models, [2]). The main point to be emphasized 
here is that there can be solutions with naked singularities when the energy in the scalar field is finite and not 
too large, and for which the effective mass of the system is positive. Thus for the simple model at hand there 
is an open set of physically acceptable scalar field data for which curvature singularities are not hidden inside 
event horizons |30U31| . This would seem to have important implications for the cosmic censorship conjecture 
[29j [38J [32] . It is worthwhile to note that, in general, naked singularities have observable consequences that 
differ from those due to black holes [37]. 

7.1 Minimal coupling to gravity 

The scalar field part of the action in curved space is 

A = \J <f*KH (i - -7= d„ (V^gg^tv) - \ g^M^j V^g d 4 x . (222) 

This gives a symmetric energy-momentum tensor Q a p for (j) upon variation of the metric. 

5A = i J J~g e aP Sg al3 d 4 x , (223) 

Q a p = 4> a <j>p (i - gT^k,) - \g aP gTW* (i - \g pa 4> P 4>a) 



9 

(224) 

It also gives the field equation for <f> upon variation of the scalar field, £ [4>] = 0, where 

5 A = - J J~g £ [(j)] Scf) d A x , (225) 
£ [0] = d a [g a ^ P V~g - 5 Q % g^M^V^g - g a %d„ (V=£<r^,) + ^y/^gg^dp WMv)] • ( 22 6) 

Since £ [<f>] is a total divergence, it easily admits a first integral for static, spherically symmetric configurations. 
Consider only those situations in the following. 

7.2 Static spherical solutions 

For such configurations the metric in generalized Schwarzschild coordinates is [35] 

(ds) 2 = e N ^ {dtf - e L ^ (dr) 2 - r 2 (dd) 2 - r 2 sin 2 9 (d V ) 2 . (227) 

Thus for static, spherically symmetric <fi, with covariantly conserved energy- momentum tensor (|224[) . Ein- 
stein's equations reduce to just a pair of coupled lst-order nonlinear equations: 



, 2 9 t * = e- L {rU - 1) + 1 , (228) 
r 2 Q r r = e - L (-rN 1 - 1) + 1 . (229) 

These are to be combined with the first integral of the 4> field equation in this situation. Defining 

77 (r) = e - L{r)/2 , w (r) = 77 (r) (r) , (230) 

that first integral becomes 

r< e -N/2 i / 4 \ 

—= W (1 + W 2 ) + -[N' + - )r l w 2 , (231) 



2 V r 
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where for asymptotically flat spacetime the constant C is given by linv-^oo r 2 4>' (r) = C . Then upon using 

t * = 0/ = = \zu 2 (1 + \zu 2 ) - riw 2 w' , (232) 
9/ = -\zu 2 (1 + \zu 2 ) - ^zu 3 (N> + ±) , (233) 

the remaining steps to follow are clear. 

First, for C ^ 0, one can eliminate N' from (I229|) and (|231[) to obtain an exact expression for N in terms 
of T), zu, and C: 

e N/2 = *C n-\rzu 3 

rzu (Azu — 2r 2 zu 3 — r 2 zu 5 + 8rrj + \2zur\ 2 + &rvo 2 r\) 

If the numerator of this last expression vanishes there is an event horizon, otherwise not. When r\ — \rzu 3 
the denominator of (|234l) is positive definite. 

Next, in addition to (|228[) one can now eliminate N from either (I229P or (|231l) to obtain two coupled 
first-order nonlinear equations for 77 and zu. These can be integrated, at least numerically. Or they can be 
used to determine analytically the large and small r behaviors, hence to see if the energy and curvature are 
finite. For example, again for asymptotically flat spacetime, it follows that 

e L ' 2 ~ i + ^- + l(6M 2 -C 2 )\ + hl (5M 2 ~2C 2 )\ + o(±) , (235) 
r— >oo r 4 ' r 2 r \t J 

e Nl2 ~ I--- Vl + lAf(C 2 -6M 2 )l + ofl) , (236) 
r— too r 2 r 12 r \ T I 



*> _~ 3 U+-+o^ 2 -? )+0 -k . (237) 



C/ M 3 „ r2 1 \ n ( 1 

-9 1 + — + o M ^ + -«! 
r — ^oo 7* \ T 2 r ' ' ^ 

for constant C and M. 

As of this writing the details of the two remaining first-order ordinary differential equations are not 
pretty, but the equations are numerically tractable. In terms of the variables defined in (|230j) . in light of 
(|234[) . Einstein's equation (|229j) becomes 

F (r, w, 77) r—zu + G (r, zu, if) r—rj = H (r, zu, if) , (238) 
dr dr 



F(r,zu,rj) = -4rj [ 2r 3 zu 6 + 3r 3 zu 8 + Wzun + Arzu 4 

+I6rr] 2 + 48zun 3 + A8rzu 2 rf + 12rzu A rf - I2r 2 zu 5 n] , (239) 
G(r,zu,ri) = 8-qzu 2 [2r 2 tu 2 + 3r 2 tu 4 - 12^ + 12rn7 3 7j + 4 ] , (240) 
H (r, zu,rf)=zu[ 8-qzu (4rtu 3 - 4j] + 2r 2 zu 2 r\ + Zr 2 zu 4 r\ + 12rzu 3 i] 2 - 12?7 3 ) 

+ (4 + 3r 2 ro 4 + 2r 2 w 2 + 12rj 2 ) (4zu ~ r 2 zu 5 - 2r 2 vu i + 8rw 2 n + 8rr] + I2mn 2 ) ] , (241) 

while Einstein's equation (I228P becomes 

I (r,zu,rj) r—zu + J (r, 07,77) r—r\ — K (r,zu,rj) , (242) 
dr dr 

I (r, tJ7, 77) = rnzu 2 , J (r, zu, rf) — —2r\ , (243) 
K (r, zu, 77) = \r 2 zu 2 (l + \va 2 ) + V 2 - 1 . (244) 



7.3 Numerical results 

As a representative example with zu > 0, (|242j) and (|238[) were integrated numerically to obtain the results 
shown in the Figure, for data initialized as zu\ r=l = 0.5 and 7?| r=1 = 1. Evidently it is true that 77(7-) 7^ 
\rzu 3 (r) for this case, so e N ( r ' does not vanish for any r > and there is no event horizon. 
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For initial values m(s)\ s=0 = 0.5 and r)(s)\ s=0 — 1.0, dip/dr — zu/rj is shown in red, e L = 1/rj 2 in green, and 
e N in blue, where r = e 3 . For comparison, Schwarzschild e L and e N are also shown as resp. green and 

blue dashed curves for the same M ~ 0.21. 

"... an exotic type of matter with which human science is entirely unfamiliar is required for such a geometry 
to exist:' — B K Tippett [34] 




e s, H (e s ) for w (s)| s=0 = 0.5 and r\ (s)| g=0 = 1.0, where r = e s . 
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However, there is a geometric singularity at r = with divergent scalar curvature: lim,.-^ r 3 / 2 i? = const. 
Since R = —Q^ 1 , and linir^o ro is finite, this divergence in R comes from the last term in (|233[) . which in turn 
comes from the second term in A, i.e. the covariant d(f>d(f>d 2 (t) in (I222p . In fact, it it not difficult to establish 
analytically for a class of solutions of the model, for which the example in the Figure is representative, the 
following limiting behavior holds. 

lim f e Lj2 l^\ = I , lim f y/Fe N/2 ) = n , lim w = p , lim U'/y/r) = pi , (245) 

where £, n, and p are constants related to the constant C in (|23 1[) : 

2C = Snp 2 /£ . (246) 

It follows that for solutions in this class, 

lim r 3/2 R = pC/n . (247) 

r— s-0 

For the example shown in the Figure: I w 1.5, n ps 0.086, p s» 3.3, C w 0.94, and pC/n pa 36. 

For the same 7?L_i = 1, further numerical results show there are also curvature singularities without 
horizons for smaller K7| , > 0, but event horizons are present for larger scalar fields (roughly when w\ , > 
2/3). A more precise and complete characterization of the data set {m\ , , 7?L =1 } for which there are naked 
singularities is in progress, but it is already evident from the preceding remarks that the set has nonzero 
measure. 

The energy contained in only the scalar field in the curved spacetime is given by 

/>oo />oo 

^Galileon - / U W dr = C S % (e S ) ds , (248) 

JQ J-oc 

H (r) = Anr 2 e L / 2 e N / 2 Q t t = 2ne 2s e L ' 2 e N l 2 vj 2 (s) (l + \vo 2 (s)) - ^e s e N ' 2 vj 2 (s) j-w (a) . (249) 

For the above numerical example, the integrand e s "H (e s ) is shown in the Figure. Evidently, -EGaiiicon is finite 
in this case. It is also clear from the Figures that the Galileon field has significant effects on the geometry 
in the vicinity of the peak of its radial energy density. There the metric coefficients are greatly distorted 
from the familiar Schwarzschild values, and as a consequence, the horizon is eliminated. 
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Other numerical examples 

Here are additional plots for r)(s)\ s _ Q — 1.0 and various initial values w(s)| s=0 . As before, d4>/dr = zu/r] is 
shown in red, e L = 1/r/ 2 in green, and e N in blue, where r — e s . For comparison, Schwarzschild e L and e N 
are also shown as resp. green and blue dashed curves for the same M, as given in the Figure labels. 




-10 -8 -6-4 -2 2 4 -10 -8 -6 -4 -2 2 4 

S S 



Initial values tu(s)\ s=0 — 0.100 and ?y(s)| s=0 = 1.00 Initial values tz(s)\ s=0 — 0.200 and v( s )\ s =o = 1-00 

corresponding to M = 0.00358 and C = 0.121 corresponding to M = 0.0191 and C = 0.283 




-10 -8 -6 -4 -2 2 4 -10 -8 -6 -4 -2 2 4 

S S 



Initial values zv(s)\ s=0 = 0.300 and r)(s)\ s=0 = 1.00 Initial values zv(s)\ s=0 = 0.400 and r)(s)\ s=0 = 1.00 

corresponding to M = 0.0546 and C = 0.484 corresponding to M = 0.117 and C = 0.710 
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Initial values zv(s)\ s=0 = 0.500 and r](s)\ s=0 = 1.00 
corresponding to M = 0.209 and C = 0.936 



Initial values zo(s)\ s=0 = 0.600 and ?7(s)| s=0 = 1.00 
corresponding to M — 0.326 and C = 1.13 
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Initial values O7(s)| s=0 = 0.700 and ??(s)| s=0 = 1.00 
corresponding to M = 0.453 and C = 1.26 



Initial values za(s)\ s=0 — 0.800 and T](s)\ s=0 = 1.00 
corresponding to M — 0.573 and C = 1.32 



For each of the last two plots, the numerical integration of the coupled galileon-GR equations has encountered 
a mathematical (as opposed to physical) singularity and terminated, resp. at r e~ 2,5 = 0.082 and 
r sa e~ 1,3 = 0.27, as is indicative of an horizon for which e N — 0. This feature persists for initial data with 
larger values of w(s)| s _ , when J?(s)| s=0 = 1. 
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Here are two more cases, just below and just above the point where horizons are formed. Again, for 
the second of these plots, the numerical integration of the coupled galileon-GR equations has encountered a 
mathematical singularity, and terminated at the point where e N ^ (blue curve) vanishes. 
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Initial values ro(s)| s=0 = 0.645 and f?(s)| s=0 = 1-00 corresponding to M — 0.383 and C = 1.199 
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Initial values zu{s)\ s=0 = 0.652 and f?(s)| s=0 = 1-00 corresponding to M = 0.392 and C = 1.209 
A useful test for an horizon is provided by the numerator of e N in (|234j) . Define the discriminant 



disc (r) = 1 — 



- ii7 (r) 
2 rjjrj 



(250) 



Should this vanish at some radius for which r\ (r) is finite, then at that radius = 0, thereby indicating 

an horizon at that radius. 

The critical case, separating solutions with naked singularities from those with event horizons, has the 
small r limiting behavior r\ (r) ~ r tn 3 (r), such that the discriminant disc ~ \ as illustrated here for 

r— > T— >0 

specific data. 
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The discriminant disc = 1 — ^rvj 3 /r] versus s = lnr for various K7| r=1 (namely, 0.4, 0.5, 0.6, 0.64, critical, 
0.66, and 0.7) with rj\ r=1 = 1. The critical initial value for the separatrix, for which disc ~ 1/2, is 

w\ r=1 = 0.65002917- ••. 

For initial data giving rise to naked singularities, disc > 1/2 (cf. the upper curves in the figure above), while 
for data leading to horizons, e N vanishes at the horizon radius, and therefore at that radius disc = (cf. 
the lower two curves in the figure). When the limiting critical behavior r/ (r) ~ r tn 3 (r) is inserted into 



the differential equations (12421) and (|238[) we find the power law behavior: 

^critical (r) ~ cr~ 3/5 , Critical (r) 



r->0 



^critical (r) ~ c 3 r 4/5 

r— f0 



„l/5 



r— >0 c 



(251) 

Moreover, critical cases are easily determined numerically for various initial data, {^(s)| s=0 , ii{ s )\s=q}^ 
thereby allowing determination of a curve that separates the open set of initial data that exhibits naked 
singularities from the set that exhibits event horizons. 

7.4 Censored and naked phases 

The situation for a portion of the initial data plane is as follows. 

H 1 1 1 1 1 



eta(r=l) 



1.2 



1.1 



0.9 



Naked singularities 




Event horizons 



o.o 



0.2 



0.4 



0.6 



varpi(r=l^° 



( w L=i ' r ?l r= i) boundary separating initial data that exhibit naked singularities from data that exhibit 
horizons. The curve is a fourth-order polynomial fit to the numerically computed critical points (dots), 
namely, 7? fit (m) = l + 0.0255538^7 - 1.34405tu 2 + 2.20589ro 3 - 0.304933zn 4 . 

This shows naked singularities for the model exist for an initial data set of non-zero measure, and are actually 
encountered for a significant portion of the initial data plane. 



40 



A similar demarcation between naked/censored solutions can be presented in terms of asymptotic r — > oo 
data instead of initial r = 1 data. With M and C denned as in (j235[) . (|236[) . and (|237[) . we find the 
following curve separating the two types of solutions. Solutions for points above the red curve have naked 
singularities, while solutions for points below that curve have event horizons. 




o-f 1 1 1 1 ■ 1 1 1 1 h 

0.0 0.2 0.4 0.6 0.8 1.0 

M 



Computed points (red circles) and an interpolating curve (solid red) separating the r — > oo asymptotic data 
for solutions with naked singularities from that for solutions with event horizons. 

By imposing the same T)(s)\ g _ initial condition for various values of t*7(s)| , the numerical data also 
shows that the corresponding C (M) has a local maximum, and hence M (C) becomes double- valued near 
that point. For example, when ?7(s)| s=0 = 1 the local maximum for C (M) is near i&(s)\ s _ « 0.78. By 
examining larger ro(s)| s=0 for the same f]{s)\ s _ , it is apparent that C (M) can also be double-valued. All 
this is evident in a parametric plot of the corresponding (C, M) points on the data plane. For example, for 
v( s )\ a =o = 1 an d various vj(s)\ s=0 £ [0.1, 1.259 92 = s/2] , we find the naked (green circle) and censored 
(black circle) data as included in the last figure, with a fitted interpolating curve (orange dashes) connecting 
the computed points. In this numerical analysis, care should be taken not to have tu(s)\ s=0 larger than 
yjl i]{ s )\s = o because otherwise this would place data initialized at r = 1 within the horizon. The horizon is 
exactly at the radius r = 1 when txr(s) | s=0 = y2 f]( s )\s=o- The S ra Y curve in the last figure is the image of 
' CJ ( s )\s=o = \/2 r)(s)\ s=0 on the (M, C) plane. Points below this gray curve can be investigated numerically 
using Schwarzshild coordinates but only if the initial data is specified for r > 1, i.e. outside the horizon. 
(Also note the portion of the initial data plane shown in the previous figure lies entirely above the curve 
T}(s)\ s=Q — i ru 3 (s) | , so all initial data points in that figure lie outside of any horizons.) 

8 Conclusions 

In conclusion, as previously emphasized by many authors it would be interesting to search for evidence 
of galileons at all distance scales, including galactic and sub-galactic, as well as cosmological. Perhaps a 
combination of trace couplings and various galileon terms, such as those in (I194[) and (|218l) extended to 
included GR effects, will ultimately lead to a realistic physical model. In particular, it is important to 
investigate the stability of galileon solutions and to consider the dynamical evolution of generic galileon and 
other matter field initial data, along the lines of [4j [5] , to determine under what physical conditions naked 
singularities are actually formed. 
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